A first-principle micromagnetic and statistical calculation of viscous remanent magnetization (VRM) in an ensemble of cubic magnetite pseudo-single domain (PSD) particles is presented. The theoretical methods developed apply to all magnetic particle sizes, from single domain to multidomain. The numerical implementation is based on a fast relaxation algorithm for finding optimal transition paths ...